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g : Abstract 

^ . In this paper, we discuss how effective environments incorporating periodic measurements can 

OO I be used to prepare a two-level system (TLS) in almost arbitrary therinal states: Concretely, we 

study a TLS coupled to a spin environment, the magnetization of which is measured periodically. In 
o ■ ensemble average these measurements cause a relaxation of the TLS into a thermal (diagonal) state. 

r^ , By adjusting the time between the measurements and the detuning of the environmental spins, 

^ I the creation of very low temperatures as well as inversion becomes possible. Our analytical results 

derived for large environments are numerically shown to be valid even for quite sinall environments, 
^ . down to only a few spins. 
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I. INTRODUCTION 

Within a quantum-therniodynaniical approach [l| a single subsystem can be shown to 
relax to a thermal state if appropriately embedded in a quantum environment: The per- 
tinent temperature is determined by the spectral density of this environment, which, for 
sufficiently large modular embeddings and for a given working point in energy space, shows 
an approximate exponential behavior. Beyond a design of this density there appears to be 
no immediate possibility of further manipulation. 

Putting aside explicitly time-dependent Hamiltonians, an interesting way to circumvent 
this limitation consists in the application of quantum measurements. Measurements play 
an outstanding role in quantum mechanics, since they provide a connection between the 
abstract theory and experimentally accessible quantities. However, in contrast to classical 
measurements, quantum measurements usually influence the measured system. A pertinent 
example for the influence of periodic measurements is the well-known Zeno effect J2|, |3|, i.e., 
the suppression of decay by fast repeated measurements. Recently, also the possibility of 
cooling a TLS coupled to an oscillator bath by applying periodic quantum non demolition 



measurements has been discussed 



Ki- 



In a recent Letter \0\ we discussed the influence of periodic measurements in a quantum- 
thermodynamical setting and pointed out how different concepts of statistical mechanics may 
arise from this "observed quantum thermodynamics". Here, we present a detailed discussion 
of the effect of periodic measurements in a generalized environment model, allowing for a 
detuning between system and environment. As we will see, this permits us to prepare 
the embedded system in almost arbitrary thermal states — even the creation of inversion 
becomes possible. 



II. MODEL 

A typical quantum-thermodynamical model consists of the system of interest and its 
environment, with which it is allowed to exchange energy via a weak coupling. The model 
we investigate here belongs to the class of systems with a modular environment. Following 8| 
we consider a system S (a TLS with energy splitting 6) described by the Pauli operator a^ 



and coupled via V (strength A) to a quantum environment B with Hamiltonian Hb, 

Htot = ^o, ® Ib + 15 ® ^B + Af . (1) 

The environment consists of n spins (e.g., a paramagnetic salt) with energy splitting 
(5^ = (5 + fi each, i.e., we allow for a detuning f2 between system and environment. We will 
show below that this detuning Q is key for preparing the TLS in almost arbitrary thermal 
states. 

Without spin-spin interaction the degeneracy of the energy level Ek = dsk with k spins up 
is given by N^ = (^). Around a given kQ, the working point in energy space, the degeneracy 
structure can be approximated by an exponential of the form 

Nk ^ Afe^""' . (2) 

For n > ko ^ 1 one finds (see Appendix [A]) 

Including a weak interaction within the environment, each energy level broadens into an 
energy band of width Ask ^ Sb- 

For the interaction between system and environment we choose 

\/ = (ff+ + a") ® 5 , (4) 

where a~^ = |1)(0| is the creation operator and a~ = |0)(1| the annihilation operator for the 
TLS and the environmental part reads 

5 = ^ ^ Ck+i,k{mk+i,nk)\nk){mk+i\ 

k nfc,mfc+i 

+ Ck,k+i{nk,mk+i)\mk+i){nk\. (5) 

To keep the model as unbiased as possible, we do not further specify the interaction but 
consider the Cij to be a set of random matrices whose entries are taken from a Gaussian 



distribution and are normalized to \Cij{a,b)\'^ = (NiNj)'^^"^. This choice of the normaliza- 
tion together with a sufficiently small A guarantees the weak coupling for arbitrary states. 
Due to Hermiticity oiV, we have Cij = Cj^^. 



Initially we assume a product state for the total system ptot(O) = /3(0) ® p%{ko), where 
the state for the environment B is given by 

Psih) = j^ '^\nko){nko\ , (6) 

i.e., only band k^ is occupied. 

Under the above conditions and for fi = the subsystem B acts as a thermalizing 
environment on the TLS: Subject to undisturbed Schrodinger evolution of the total system 



(i.e., without measurements), the TLS would relax into a state wit 



Hk 



1 inverse temperature 



io|. 



T"^ = /3(fco) given by (|3]) independently of its initial state /3(0) 

In the following, we want to study the effect of external observation (by means of periodic 
measurements) for this quantum-thermodynamical resonant case as well as for the case of a 
detuned environment (i.e., 6b t^ S). 

III. EFFECTIVE ENVIRONMENT 

Without measurements, the environmental spins have to be in resonance with the TLS, 
i.e., the detuning has to vanish {Q = 0) in order to allow for a thermalization of the TLS. 
However, this no longer needs to be the case if we include periodic measurements. We 
thus analyze the dynamics of the system with arbitrary detuning Q disturbed by periodic 
measurements of the environmental magnetization, i.e., the system first evolves under pure 
Schrodinger dynamics for some time At. Then a measurement of the environmental mag- 
netization (or rather its energy) is executed, which means we perform a projection on one 
of the energy bands ki. The pertinent projection operator reads 

Ai = 15 ® Pk, (7) 

with 

Ai = ^|^fei)(nfcj , (8) 

where we sum over all levels of the measured energy band ki. This leads to a state of the 
TLS given by 

P(l) = fz f, (9) 

Trs|Pfc,p^(l)| 
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which can be rewritten as 



p(i) = p'(i) + — y^ ^, (10) 



where p' denotes the state of the TLS before the measurement and the second term is the 



so-called co-jump [ll| caused by the correlations Csb '■= Ptot ~ p' ® p's between system and 
environment straight before the measurement. Since the measurement is incomplete (mea- 
surement of an energy band and not of a single level), there may still be some correlations 
left after the measurement. However, the effect of these correlations is negligible in the 
case of weak interaction and small band width. That is, the state of the total system after 
the measurement can be approximated by a product state. For simplicity, we introduce a 
further approximation: Since we are not interested in details like the occupation probabili- 
ties of the single levels within the energy bands of the environment, we apply some kind of 
coarse graining by replacing the typically rather complicated state of the environment after 
the measurement by a simple one of the form (jH]). This can be done, provided the environ- 
ment is sufficiently large. Thus, the new state of the total system after the measurement is 
approximated by 

Ptot(l)^p(l)®P°B(fcl), (11) 

i.e., a form similar to the initial total state but with a new state p(l) for the system and 
possibly a different energy band being occupied for the environment (if a different energy 
band has been measured). This procedure can now be repeated periodically. In order to 
analyze what happens with the TLS due to these measurements, we first have to calculate the 
state of the TLS after a single measurement dH]). This can be done by means of perturbation 
theory (see Appendix |B]), which leads to three possible states of the TLS corresponding to 
the three possible measurement results: 

In the case of measuring the same energy as in the previous measurement, we obtain the 
occupation probability of the ground state up to second order in the interaction strength A: 

Poo,o(j) ^ Poo(j - 1) (1 - 4AVii(j - l)(e-^^-/2 - e^'-/') (Ci - C2)) (12) 

with 



sm 



2 (n 



(fA() 



C, := ]^, ' , (13) 

^' (2(5 + ^)2 ^ ' 



In case of measuring one energy band higher than before, we get 



^■^_ Pii(j-l)Ci /,KN 

PiiU -1 Ci + PooU - 1 C2 



whereas measuring one energy band lower leads to 



Poo,- J ~ J-. 7T-— J-. TT-r ■ (16) 

Pii u - 1 C2 + P00 u - 1 Cl 



For the off-diagonal elements we get accordingly 



Pio,o(j) ~ Pio(j - 1) 



l + \H (e-^^-/Voo(j - 1) + e^^-/Vii(j - 1)) 4Ci 



1 - e^^^* + ifiAt 1 - 6^(2-5+?^)^* + i(25 + fi)At\ 
W- ^ {25 + nf J 



(17) 



, 1 + e^'^^' - 2&^^' cos((5 + fi) At) 
4(2M2 + 12^) (pii(j - l)Ci + PoolJ - 1)C2) 
l + e2^^-"*-2e^^-"*cos((^ + l])At) 

Pio,- U -Poi J - ^> Afnxn^n2\f — T — Tv^^ T — TvM ' (^^) 

4(2M2 + i2^) (pii(j - 1)C2 + PoolJ - l)Ci) 

As one can see, measuring a different band causes significant changes in the occupation 
probabilities of the TLS: For example, if poo was close to 1 at the previous measurement it 
will be almost after such a measurement, and vice versa. Thus, one cannot expect that a 
single TLS will relax to a stable thermal attractor state as is the case without measurements 
in the quantum thermodynamical model. Indeed, after sufficiently many measurements the 
TLS will always be in its ground or its excited state, i.e., passing through some kind of 
quasi-classical trajectory. 

Instead, one may ask what happens to the TLS in ensemble average. Therefore, we 
first have to weight each possible outcome with the probability to obtain the corresponding 
measurement result. These probabilities are given by (]B14p : 

which equals the occupation probabilities of the bands kj = {/cj-i, %-i + 1, kj-i — 1} before 
the measurement. 

Thus, for a given energy band kj_i determined at the preceding measurement j — 1, the 
probabilities for measuring one band higher, one band lower, or the same band, respectively, 



read 

p+(j) ^ 4AV^-/2 (pn(j - l)Ci + Poo(j - 1)C2) , (21) 

pM) ^ 4A2e-^^-/2 (poo(j - l)Ci + PiiU - 1)C2) , (22) 

Po(j)~l-P+(j)-P-(j)- (23) 
This yields the ensemble average (indicated by an overbar) for the diagonal elements, 

Poo(i) = Po{j)poo,o{j) + P+{J)pooAJ) + P-iJ)poo,-U) 

^ PooU - 1) + 4A2 [ {e^'-/'pnU - 1) - e-^^-/Voo(j - 1)) Ci 

+ (e-^^-/Vii(j - 1) - e^^-/Voo(j - 1)) C2] . (24) 

We can now iterate this result to get the ensemble-averaged density matrix as a function 
of the number of measurements, i.e., the ensemble average after measurement j is used as 
initial state for measurement j ' + 1 and so forth. The calculation for the off-diagonal elements 
can be found in Appendix O It turns out that the off-diagonal elements vanish for almost 
any choice of the parameters after sufficiently many measurements (an interesting exception 
is discussed below), i.e., the attractor state of the TLS will be a thermal state. We note in 
passing that the ensemble average coincides with the infinite-time average (ergodicity). 

IV. EFFECTIVE RELAXATION CONTROLLED BY At AND Q 

To derive the analytical expression for 'PooiJ)^ "we first rewrite (|24l) as 

Poo(j + 1) - Poo(j) = -RpooU) + d (25) 

with 

R := SA^ cos\i{^6b/2) (Ci + C2) , (26) 

rf:=4A2(e^^^/2Ci + e-^'^/%) . (27) 

The left-hand side can be approximated by a derivative with respect to j, as long as pgo 
changes only slowly with j, which is guaranteed by the weak coupling (A ^ 1). Thus, we 
can approximate ( l25l) by a differential equation of the form 

^^ = -RPooU) + d, (28) 



with the solution 



Poo(j)=(Poo(0)-|)e-«^ + |, (29) 



This means that the ensemble of TLSs exponentially approaches an attractor state with 
increasing number j of measurements. 

This relaxation has some interesting properties: First, it is noteworthy that due to the 
measurements a relaxation becomes possible not only in the resonant case (as without mea- 
surements) but also for 6b 7^ S. However, the attractor state reached as well as the relaxation 
constant depend on the time At between the measurements and on the detuning Q. 

Considering the relaxation constant 

, _ /sin2(Mt) sin2((5+^)At) , , , 

R = 8A^ coshW6s/2) { -^^ + ^^^^ ^ ^^, ^ ] , (30) 



one finds 



lim i? = , (31) 



which means that the relaxation is slowed down for very rapidly repeated measurements. 
Such a suppression of decay due to fast, periodic measurements is well known as the so- 
called quantum Zeno effect J2|, yj. Furthermore, R also vanishes for At = n^ and Q = ^^, 
n = 1,2,..., m = 1,2,.... For this special choice of the parameters (characterizing the 
exceptions we mentioned before), the absolute value of the off-diagonal elements also stays 
constant, as discussed in Appendix O Thus, it is possible to freeze the initial state of 
the TLS and suppress decoherence by periodic measurements carried out even with finite 
frequency. 

For any other choice of the parameters, however, a thermal attractor state is obtained, 
which also shows remarkable features. In general, it is given by (cf. (1291) ) 

d 



— attr 

Poo -^ 



QfiSB/2^^ + q-PSb/2(^^ 



2 



2cosh(/35B/2)(Ci + C2) 
Let us first study the resonant case fi = 0. Because of 



(32) 



sin^ (§At) At^ , , 

hm Ci = lim ^^^ = -— , 33 

n^o n^o Q? A 

lim C - lim ^^^M(^ + f)^^) _ '^' i^^t) (u) 



we get 

^°° ^ ' ~ 2 cosh(/3V2) {sm\6At) + 6' At') ' ^ ^ 

The lowest possible temperature for the TLS in this resonant case is obtained for At = ^, 
n = 1, 2, 3, . . . and is given by Tmin,res = 1//3- This is just the temperature the TLS would 

n 

also get according to quantum thermodynamics without measurements (cf. [8]). However, 
any other choice of At will lead to a higher temperature. In particular, for rapidly repeated 
measurements we obtain 

iim#' = i. (36) 

i.e., the TLS is heated up to very high temperatures {T^s — ?■ C)o). 

Finally we want to focus on the effect of a detuning Q ^ between system and environ- 
ment: Now, the attractor state fl32l) reads 



-attr ^ ^" "^ {2S+nr /oyx 

2cosH/3Sb/2) (^— ^ + (^^W^)^ 

and has two new interesting properties compared to the resonant case discussed before. The 
lowest possible temperature here is reached for At = m:/{5 + fi/2), n = 1,2,3,..., and is 
given by 

Ttnm = T-13- (38) 

Ob P 
This temperature depends on the ratio of the energy splittings of the TLS and environ- 
mental spins, that is, the off-resonance may be used for cooling the TLS down to very low 
temperatures by choosing proper values for At and Q. 

The other remarkable feature is the limit of the lowest possible occupation probability of 
the ground state, i.e., the highest temperature. As in the resonant case, rapidly repeated 
measurements lead to Tes — )■ oo. However, in the off-resonant case, it is even possible to 
further increase the energy of the TLS, i.e., to create inversion. The maximum inversion is 
reached for At = 2mr/\Q\, n = 1,2,3, .. ., and is given by 

p-/35s/2 

^^iPor(^t)) = ,,s./2 + ,-,s./2 - (39) 

which corresponds to the negative temperature [12] 

T^ax = -^^. (40) 

Ob P 




FIG. 1. Attractor for the occupation probability of the ground state Poo*"^ ^^ ^ function of the 
parameters At and (7. States with low temperature are white and those with high (negative) 
temperature are black. Parameters used: 5 = lu, /3 = 0.75 n~^. 

Note, that this inversion originates from the influence of the periodic measurements and not 
from a negative temperature state of the spin environment (the environmental /3 is positive). 
By varying the time At between the measurements, thermal states of arbitrary tempera- 
ture within these limits can be prepared. Possible values of inverse temperatures thus cover 
the range between —{5b/5)I3 and +{5b/5)[3, symmetrically around zero. (In the resonant 
case these values are between zero and /3 only.) Fig. [1] shows the attractor for the occupation 
probability of the ground state as a function of time At and detuning Vt. 

V. NUMERICAL VERIFICATION FOR SMALL ENVIRONMENTS 



The discussed behavior of the TLS under periodic measurements of the environment can 
be obtained not only for large environments but even for quite small ones consisting of a 
few spins only. For such small environments the attractor state for the TLS depends not 
only on the initial state of the environment (working point) but also on the initial state of 
the system, since the environmental degeneracy structure can no longer be approximated 
by an exponential behavior. This means that the system no longer shows thermodynamic 
behavior, which requires independence of its initial state. Nevertheless, given the initial 
state, we can predict the attractor state using our analytical results. 

To demonstrate this, we consider a TLS with splitting 5 coupled via a random [a^ ® cFx)- 
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FIG. 2. (Color online) Relaxation of the TLS given by the ensemble-averaged pgo under periodic 
measurements of its environment (seven spins with initially two spins up) . We obtain a relaxation 
to the thermal attractor state with 'p^^^ = 3/4 (dashed line). Chosen parameters: 5 = lu, Q = 0, 

At = tt/6. 



interaction to a seven-spin environment (splitting of the spins 6b), with initially two spins 
up. Let us first consider the resonant case 6 = 6b = lu with At = 7i/6 and the TLS initially 
being in its ground state. In this case only the initially occupied environmental "energy 
band" E2 = 26b and the band beneath (£"1 = 6b) are involved in the dynamics (the ground 
state energy is set to Eq = 0). This means that only the ratio of the degeneracies of these 
bands determines the final attractor state. Thus, because A^^i = d) = 7 and N2 = (2) = 21, 
we get an effective inverse temperature /3eff = ln{N2/Ni)6j^^ = ln(3)5^^ and therefore expect 
Pqq*'' = 3/4. Indeed, this final state is obtained in our numerical exact simulation shown in 
Fig. El 



In the off-resonant case, it is possible to obtain inversion as discussed in the last section. 
To illustrate this, we choose Q = 0.7 and At = 2tt/Q. Again, initially we start with two spins 
up. Since now only the second and third bands are involved, we obtain (3es = ln(35/21)5^^, 
which yields according to (139|) Pqq*'' = 3/8. Again, this value is verified by the simulation, as 
shown in Fig. [3l 
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FIG. 3. (Color online) Relaxation of the TLS to a state of negative temperature due to the 
periodic measurements of its environment (seven spins with initially two spins up). The final 
state is determined by /Jqq*'' = 3/8 (dashed green line). Chosen parameters: 5 = lu, i^ = 0.7 u, 

At = 2TT/n. 

VI. CONCLUSION 



In this paper, we have studied the influence of an effective environment on a coupled two- 
level system. We have shown analytically that periodic measurements of the environmental 
magnetization (energy) may lead to a freezing of the initial state, i.e., suppress decoherence, 
or may (in ensemble average) yield a thermal state for the TLS, depending on the choice 
of the detuning and the time between the measurements. The temperature of the attractor 
state reached is controlled by these two parameters, which allows for the preparation of 
thermal states with almost arbitrary temperature. It is even possible to create an inversion 
for the TLS. The key mechanism, which underlies this cooling or creation of inversion, is the 
periodic destruction of the correlations between the system and its detuned environment. 
The resulting periodic reset of the interaction energy allows for a successive change of energy 
of the system and environment, even if out of resonance. It is to be expected that the 
deviation Tes from T would give way to non-thermal attractor states in the case of larger 
subsystems S, i.e., beyond the two- level case. [l3| 

As demonstrated in our numerical simulations, this intriguing behavior can be observed 
even for quite small environments consisting of a few spins only. In this case, the attractor 
exhibits dependence on the initial state of the system S. 
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In a recent paper [IJ] a nuclear magnetic resonance setup has been used for studying the 
influence of periodic quantum nondemolition measurements in a related model. We think 
that similar experiments may also be suitable to experimentally test the results we have 
derived here. 

Appendix A: Thermalizing environment 

In this appendix, we show how the degeneracy structure of a system consisting of many 
spins can be approximated by an exponentially increasing degeneracy around a certain 
point in energy space, which is a precondition for this system to act as a thermalizing 
environment |ll|]. Thus, we want to approximate 

"f] ^ Afe"'^" (Al) 

around some working point ko (we set the ground state energy to 0). To do so, we first 
replace the discrete binomial coefficient by a continuous function, which can be done using 
n! = r(n + 1), with T{x) being the Gamma function. Thus, the binomial coefficient can be 
written as 

"^ = r(" + ^) f A2) 

Now, we can apply a Taylor expansion around ko, which reads 

:) = nn-kHmh^i) <^ ^ «'<" - *° ^ ^' - *<*° ^ ^>"* - *°» ^ '^ «* - *°*'' 

(A3) 
where ■ip{x) := ^ln(r(a;)) is the so-called Digamma function [15]. Comparison of this with 
the Taylor expansion of Me^^'^^ yields in zeroth order 



ko\{n- /co)! 
and in first order 



j^e^^Bko (A4) 



ko\{n-ko)\ 
Hence, with fjA5l)/flA4l) we get 



""■ {iP{n - fco + 1) - ^(fco + 1)) = Mp^Be^^^''' . (A5) 



P5b = ij{n -ko + l)- tPiko + 1) (A6) 
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For large x the Digamma function can be approximated by ip{x) ~ ln(x), which finally leads 
to equation Q, i.e., 

This means, that such an environment consisting of many spins with total energy Ssko can 
act as a thermalizing environment on a coupled TLS with 6 = 6b, enforcing a thermal state 
of temperature 1//3 on the TLS according to (]A7p . 



Appendix B: Calculation of the short-time dynamics 



The von Neumann equation for the total system in the interaction picture (we set ^ = 1) 
reads 



d_ 

dt' 



Ptot(t) = i Ptotit),XV{t) 



with 



and 



V{t) = d+ B{t) + o- B\t) 



B{t) = e^HBt-^^^iHBt^iSt 

k nk,mk+i 



fBll 



(B2) 



(B3) 



where uj{mk+i,nk) is the energy difference of level m in band k + 1 and level n in band k 
minus 6, i.e., always much smaller than S for small band width. 

Since our goal is to study the infiuence of the periodic measurements, we denote the state 
of the system by the discrete number of the measurement instead of the continuous time, 
i.e., p{j) = pit = jAt) and so forth. According to (JBip . the density operator of the total 
system up to second order in interaction strength can then be written as 



P'totU) ^ PtotU - 1) + iA / Ptot(j - 1), Vif) 



At 



dt' 



-A^ 



At rt' 



pUj-i),vit") ,vit') 



df'dt' 



(B4) 



'0 ^0 

One can easily show that the first-order term does not contribute to the state of the TLS 
since 



Trs{pB(j - l)B{t')] = Ttb{pb{j - l)^t(i')} = 0- 
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The commutator appearing in the second-order term reads 

■ptot(j-l),nt")],m = 

Ptot(j - l)a^Bit")a-B^t') + ptot(j - l)a-B\t")a+Bit') 
-a+Bit")pUj - l)o^B{t') - a^B{t")Ptot{j - l)a-B\t') 
-a-B\t")pU3 - l)a^B{t') - ^-5t(t")ptot(j - l)a-B\t') 
-a+B{t')ptot{j - l)a+B{t") - a-B\t')pU3 " l)^+5(t") 
-a+5(t')ptot(j - l)a-B\t") - a-B\t')p,,,{3 - l)a-B\t") 

+a-B\t')a+Bit")pUj - 1) + a+5(t')a-5t(t")ptot(j - 1) 
=: c. 

This leads to the 00 component of the numerator Ttb\ Pk^p'totiJ) r in (E] 



(B5) 



/At pt' 
J PooU-l)Ut\t") 

-PiiU-iMt',t"Wdt' 



(B6) 



with 



Ut',t") := TiB{Pk,PBU - l)B\t")B{t') + P,^B\t')B{t")pB{j - 1)} 






Nk 



■' mfc^,nfc^,_i 



X cos((ll + w(mfc nfc i))(t' - t")) 



AT., 



7 ^ |C'fej,fcj+i("^fcj,'"'fcj+i 



■J ' '^j 

2 



™fc^:"fcj+l 



ei(t',t" 



X cos((2(5 + fi + wK^.+i, mfcj)(t' - t")) 
TrB<! Pk,B\f)pBU - l)5(t') + A 5Ht')PB(j - l)B{t" 



> , |Cfc^,fc^._i(mfcj.,nfc^_i)| 



iVfc,_i 



^ - mk.,nk.--i 



25fej_i-lA;j 



iVfc,+i 



X cos((n + a;(mfcj., rafc^,_i))(t' - t")) 

/ , |C'A:j,fcj+l("^fcj,''^fcj+l)| 



J ■ - mfe^,nfc^.+i 



X cos((25 + fi + a;(nfc^.+i, mfcj)(t' - t")) , (B7) 
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where Sab denotes the Kronecker delta, ^o and ^i both consist of two double sums over 
many oscillating terms with normally distributed amplitudes and approximately uniformly 
distributed frequencies. For sufficiently small width of the bands (or rather short enough 
times) these sums can be approximated by 

J2 \Ca,bima, nb)\^ cos((a + a;(m„ n^W - t")) ^ Ar,iVfe|Ca,6(ma, n^) ^ cos(a(t' - 1")) (B8) 

with a being Vt oi 25 + fi, respectively. Using \C~{a~hj)^ = {NiNjY^^'^ = e~'3'^s(i+i)/2^ ^^ 
obtain 

6(t', t") ^ 24,_,fc, (e-^'-/' cos{Q{t' - t")) 
+e^^^/2cos((25 + fi)(t'-t"))), 
6(t', t") ^ 24,_,+ifc,e^^-/2 cos(fi(t' - t")) 



+2(5fe,_,_i,,e-^'^«/^ cos((25 + fi)(t' - t")) . 



(B9) 



With this, the integrals in flB6P can be easily calculated, leading to 



Trs |Pfc,Ptot } ~ PooiJ - l)4,.ifc,. 



+4A2 - PooiJ - 1)5.,_,., (e-^^-/^Ci + e^'^/\,) 
+Pii(j-1) 

x(e^^-/\i4,_,+i.,+e-^^-/^C24,_,-i.,) 



where we have introduced the abbreviations 



(BIO) 



Ci:= 
(2: = 



sm 



2 rn 



(fAt) 



fi2 ' 

sin^((^ + f)At) 
{26 + Qy 

16 



(Bll) 
(B12) 



The off-diagonal element Tr^ I PkPtot \ can be calculated in the same way, leading to 



Trij<^ PkJ'totU) \ - PwU - l)4.,_ifc 



10 



J -J- "-J 

At rt' 



-A- 



I I Pio(j - l)TrB{A,Pij(j - l)5t(t")5(t') 

-poi(j - l))TTB{Pk,Bit")pBij - l)B{t') 

+PkM')pB{3 - l)B{t") \dt"dt' 



PwU - l)'^fe,-ife, - -^^ 



4,_,., (e^^-/^ + e-^'-/') 



1 - e'^^* + i^At 1 - 6^(2^+^)^* + i{25 + n)At\ 
^ ' f]2 + {26 + Qy ) 

xpioij - 1) 

^ 26n + fi2 



xpoi(j-l) 



(B13) 



To get the density matrix of the TLS, we finally have to calculate the denominator in 
which yields 



T:rB{Pk,p'BiJ)} ^ S, 



^j - 1 ^j 



X' 



At ft' 



^0 



Trs 4Tr5{c} dt"dt 



^4,_,fc, -4A2 



Pn(j-l)(e^^-/^Ci + e-^^-%) 



+Pooij-l){{e-''^^'Ci + e^'^/%) 



Sk 



_p/3<5s/2 



(pii(j - i)Ci + Poo(j - 1)C2) 4,_i 



+lfc,' 



_„-/35s/2 



(P00(j - l)Cl + Pll(j - 1)C2) 4,_i-lfc,- 



(B14) 



Equations ( IBIOI) . ( JB13I1 . and (1B14I) thus lead to three different values for poo(j) and pio(j) 
after the measurement (equations (fT2|) - (TT9|) ). depending on whether a lower, the same, or a 
higher energy has been measured compared to the previous measurement. 
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Appendix C: Dynamics of the off-diagonal elements 

In this appendix we derive the dynamics of the off-diagonal elements of the TLS state 
due to the periodic measurements. Using (TT71) - (IT9|) as well as (12T]) -( 123|) . we find for the 
ensemble average 

Pio(j) - Pio(j - 1) + 2A' cosh{(36B/2) 
X -Pio(j-l) 

1 - e*^^* + iQAt 1 - eH25+f^)At + [(^26 + Q)At\ 

IT^ ^ {26 + ny ) 

+^°i^-^ ~ ^^ 26Q + Q^ J ' ^^^^ 

i.e., the dynamics of the off-diagonal elements can be treated independently of the diagonal 
elements. We first rewrite ( JCip as 

Pio(j) - Pio(j - 1) = Pio(j - l)(ci + ic2) + Poi(J - 1)(C3 + ic4) (C2) 



with 

c, := -2X^ cosh(/35,/2) ^ ^^^ + ^^^^^^} (C3) 

C2 := -2A2 cosh(/3(5B/2) 

/f]At - sin(l]At) (2(5 + f])At-sin((25 + fi)At)\ 

"^ ^ n^ + {26 + ny ) ' ^^^^ 

.x2 ..«. ,^,cos((5At)cos(((5 + fi)At)-cos2(5At) ,^, 

C3 := -4A^cosh(/35,/2)^ ^^ , (C5) 

Due to the weak coupling, we can approximate the left-hand side of ( JC2I) by a derivative, 
which leads us to two coupled differential equations for 9^(j) := Re(pio(j)) and J(j) := 
Im(pio(i)): 

d^(^) ^ (c, + C3)9^(j-) + (C4 - C2) J(j) , (C7) 



dm 



Ci - Cs)3{j) + (C2 + C4)^(j) . (C8) 



dj 
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The solution of this system of differential equations reads 

1 



Hj) 



:?(j) 



^{ci-7)j 



(c3(e2^^-l)+7(e'^' + l))5^(0) 



27 

+ (C4 - C2) (e^^^' - 1) 3{0) 

^e^'^^-^)^' [ (C2 + C4) (e^^^' - 1) 9^(0) 
+ (-C3 (e'^^' - 1) + 7 (e'^^' + 1)) 3f(0) 



(C9) 



(CIO) 



where we introduced the abbreviation 7 := y— cf+cf+cf. This finally yields the absolute 
value of pio(j) 



3CIJ 



7 



((c3sinh(7J) + 7cosh(7J)) D^(0) + {04 - C2) sinh(7J)Df(0))' 



+ ((C2 + C4) sinh(7j)9^(0) + (-C3 sinh(7j) + 7cosh(7j)) 3(0))' 



Considering the special case At = rij and fi = ^^, n = 1, 2, . . ., m = 1, 2, 



1/2 



that C2 = -A^ cosh{f36B/2) 



n^(2m+n)7r 
m{m+n)(5^ 



and Ci^3^4 = 0. This directly leads us to 



(Cll) 
we find 

(C12) 



i.e., a slow oscillation of the off-diagonals. 

Thus, the absolute value |pio(j)| stays constant for this special case, for which also the 
relaxation constant R fl30l) vanishes. For any other choice of the parameters, we obtain 
from (JClip limj^oo |Pio(i)l = 0; i-^-, the TLS relaxes into a thermal state. 
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